Faster quantum walk algorithm for the two dimensional spatial search 
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We consider the problem of finding a desired item out of A'^ items arranged on the sites of a 
two-dimensional lattice of size \/N x \/iV. The previous quantum walk based algorithms take 
0{y/N\ogN) steps to solve this problem, and it is an open question whether the performance can 
be improved. We present a new algorithm which solves the problem in 0{\/ N log A'') steps, thus 
giving an 0(-^log N) improvement over the known algorithms. The improvement is achieved by 
controlling the quantum walk on the lattice using an ancilla qubit. 
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I. INTRODUCTION 



Suppose we have N items arranged on a two- 
dimensional lattice of size \/N x ^/n . Let the sites 
be labeled by their x and y coordinates as \x^y) for 
x,y & {0, . . . , ViV — 1}. In the quantum scenario, the 
coordinates label the basis states of an N dimensional 
Hilbcrt space. Let f{x, y) be a binary function which is 
1 if the item placed on \x,y) site satisfies certain proper- 
ties (i.e. is a marked item m), else it is 0. We assume that 
there is a unique marked item and let \m) = \x, y) f(x,y)=i 
denote the corresponding site or basis state. The two di- 
mensional spatial search problem is to find \m) using min- 
imum time steps, with the constraints that in one time 
step we can either examine the current site (i.e. compute 
/(x, y) for the current site using an oracle query) or move 
to a neighboring site. 

The straightforward application of Grover's search al- 
gorithm [l| cannot be used to solve the problem faster 
than classical search as pointed out by Benioff Al- 
though it can find |m) using 0{^/N) oracle queries, be- 
tween successive queries, it needs to perform a reflection 
about a superposition of all sites. This reflection takes 
0{-\/N) time steps, as in one time step we can only move 
to a neighboring site and we have to move across ^/N 
sites in each direction of the lattice to perform a reflec- 
tion. Note that in the standard search problem, there is 
no restriction on the movement on lattice, and hence this 
reflection is not a hurdle. But for 2c? spatial search, the 
total complexity becomes 0{-\fN x ^/N) — 0{N) time 
steps, no better than brute-force searching. 

Aaronson and Ambainis have shown that using a clev- 
erly designed recursion of the quantum search algo- 
rithm, the 2d spatial search problem can be solved in 
©(VlVlog^ N) time steps [3]. A better alternative is pro- 
vided by the quantum walk search algorithms. They have 
been constructed for spatial search in any number of di- 
mensions (see, for example, Refs. 0, H, S 0] ) ■ For the 2d 
spatial search problem, the discrete time algorithm by 
Ambainis, Kempe, Rivosh (AKR) and the continuous 



time algorithm by Childs and Goldstone (CO) can do 
the job in 0{^/N log N) time steps. It is an open question 
whether the algorithms can be further improved, par- 
ticularly whether the lower bound of n{\/N) Q can be 
achieved. Here, we give a positive answer to this question 
by presenting an improved algorithm that can solve the 
two-dimensional spatial search problem in 0{y/N log N) 
time steps, thus giving an 0{^/logN) improvement over 
the best known algorithms. 

We present our results in the context of AKR's dis- 
crete time quantum walk algorithm, but the same can be 
applied to the continuous time quantum walk algorithm 
of CG. These quantum walk algorithms start with a uni- 
form superposition of all sites and achieve a particular 
state, denoted by la"*") here, in 0{y/N log N) time steps. 
The overlap of |a+) with |m) state is 9(l/-\/IoglV), so 
that 0{\/\ogN) rounds of quantum amplitude amplifi- 
cation 9] can be used to get the \m) state with con- 
stant probability. Hence, the total complexity of the al- 



gorithms is 0{^/N log N X Vlog N) = 0{VNlogN). In 
the case of AKR's algorithm, the quantum walk search 
is analysed by reducing it to an instance of the abstract 
search algorithm, which is a generalization of Grover's 
search algorithm. 

We modify the quantum walk algorithms in a par- 



ticular way so that the 



state, obtained after 
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0{^/N log N) walk steps on the uniform superposition, 
has a significant overlap with the |to) state. Hence no 
rounds of quantum amplitude amplification are required 
by the new algorithm, and the time complexity remains 
0{^/NlogN). As we show, this improvement is possible 
by controlling the quantum walk on the lattice in a clever 
way using an ancilla qubit. Our algorithm applies to any 
instance of abstract search algorithm, so it can also be 
used for improving the spatial search in higher dimen- 
sions. But there the improvement is only by a constant 
factor. 

The paper is organized as follows: In section II, we re- 
view the abstract search algorithm, presented by AKR, 
with the example of two-dimensional spatial search. Al- 
though our analysis follows AKR's paper, we use different 
notation for convenience. In section HI, we present the 
controlled quantum walk algorithm. We conclude the 
paper with some discussions in section IV. In the Ap- 
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pendix, we present analysis of the abstract search algo- 
rithnijwhich closely follows AKR's analysis (see section 
7 of 0) and uses the results presented there. The dif- 
ference is a minor modification which is required for the 
new algorithm. 



II. BACKGROUND 

Grover's search algorithm starts with an initial state 
\s), normally chosen to be uniform superposition of all 
the basis states. The algorithm drives it to the target 
state \t) by successively applying the reflection operators, 
Rt = 2\t){t\ — In and Rg = 2\s){s\ — In, where In is the 
A^-dimensional identity operator. The \t) (|s)) state is an 
eigenstate of reflection operator Rt (Rs) with eigenvalue 
1, and all the states orthogonal to \t) (|s)) have eigenvalue 
— 1. It has been shown that applying the operator Ug — 
Rs ■ Rt on |s) rotates it in the two-dimensional subspace 
spanned by \s) and \t), and after 0(l/|(t|s)|) iterations 
of Ug we come very close to the \t) state. 

The abstract quantum search algorithm is a general- 
ization of Grover's search algorithm, where the operator 
Rt remains the same but Rs gets replaced by a more gen- 
eral operator U. The |s) state is still required to be an 
eigenstate of U with eigenvalue 1, but the states orthog- 
onal to \s) need not be its eigenstates with eigenvalue —1 
as in the case of Rg- Also, U is required to be a real op- 
erator (not necessarily a reflection operator) and not to 
have any other eigenstate with eigenvalue 1 apart from 
\s). The abstract search algorithm iterates the operator 
Ua = U ■ Rt to get to the target state \t). 

To quantify the number of iterations of Ua needed to 
get to the \t) state, we note that since J7 is a real unitary 
matrix, its non-±l eigenvalues come in pairs of complex 
conjugate numbers e^**. The eigenstate corresponding to 
eigenvalue 1 {6 = 0) is |s), also denoted as |$o) here. Let 
the eigenstates corresponding to eigenvalue —1 {9 = it) 
be denoted by \^k), fc = 1,...M. Let \^f) denote all 
other eigenvectors with non-±l eigenvalues e^^^^ . Then 
|$+) = 1$-)* as U is real. Let af = {^f\t), ao = ($o|t) 
and Ofc = {^k\t) be the expansion coefhcients of \t) in the 
eigenbasis of U. Since \t) is a real vector, = {aj)*, 
and upto a global phase, \^f) can be chosen such that 
O'j^ — 0.J — Similarly upto a global phase |<i>o) and 
|$fe) can be chosen such that aq, Uk are real. Thus 



\t) = ao|$o) +J2a,{\'f+) + |cf>T)) + ^a,|$, 



(1) 



To analyse the iteration of operator Ua — U ■ Rt on 
|$o): its eigenspectrum was determined by AKR. Though 
they have not explicitly considered the possibility when 
U has an eigenspace with eigenvalue —1, their analysis 
can be easily generalized to that case. As it is crucial 
for the new algorithm, we have done this analysis in Ap- 
pendix for completeness. For particular cases (spatial 
search is one of them) , only two eigenvectors | ± a) of Ua 



with the eigenvalues e^*" are important, with the start- 
ing state |<i>o) almost completely spanned by them. Here 
a depends upon the eigenspectrum of U as 
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(2) 



where Ak = y Efcli '^l is the projection of \t) on the — 1- 
eigenspace. As shown by AKR, |$o) is close to \a^) = 



V2 



a 



a)). Quantitatively 



l(<i>o| 



(1 — cosOjY 



-e 



AW 



(3) 

After iterating the operator Ua for T — [7r/2a] times 
on |q;~), we come very close to the state |a"*") = 
-i(e"/2|a) - e-*"/2| _ a))/V2 = (|a) + | - a))/^/2. As 
shown by AKR, the quantity |(i|Q;''')| depends upon the 
eigenspectrum of U as 



a' 



e 



1 



(4) 



Consequently, any operator U can be used in place of 
Rs for quantum search if it satisfies the conditions for 
abstract search algorithm, i.e. it is a real operator with 
the initial state \s) as its unique eigenstate of eigenvalue 
1. Sometimes this flexibility is very useful. In the case of 
spatial search, there is a restriction that in one time step, 
we can move only to neighboring lattice sites. In this 
case, U can be chosen such that it can be implemented 
in only one time step, whereas Rg takes 0(v^) steps. 
For any U , we need to find its eigenspectrum, and the 
expansion coefficients of the target state in its eigenbasis, 
in order to analyse the algorithm using Eqns. ([l]|4]). 



A. Two-dimensional spatial search 

We illustrate the abstract search algorithm with the 
specific example of two dimensional spatial search. 
AKR's algorithm attaches a 4-dimensional coin space Tic 
to the Hilbert space TLn associated with N lattice sites, 
and works in the joint Hilbert space Tij — Hc^Hn- The 
four basis states d = 0, 1, 2, 3 of He represent the four pos- 
sible directions of movements on a two-dimensional lat- 
tice, i.e. I ^),| ^),| t),U). Let = i^Jd) be their 



uniform superposition and let \un) — Ek y y))/VjV 



be the uniform superposition of all lattice sites. The ini- 
tial state |$o) of AKR's algorithm is \^o)akr = \uc}\un) 
which can be prepared in 2y/N time steps. (For prepar- 
ing \un), the idea is to start with a site |0, 0), first spread 
the amplitude along x axis in a/TV steps, and then repeat 
the process for y-axis in another y/N steps 
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The algorithm then iteratively applies the operator 
Uw 



= W ■ Ruc.m on |$o)- The operator R^c 



-Rii 



L 



iN - 



-2\uc,m) {uc, m\ is the negative of the re- 



flection about the |uc)|"T') state. It can be implemented 
in one time step by examining the lattice sites (in quan- 
tum superposition) using an oracle, and then applying 
Rue = —Rue = li — 2|uc)(mc| iff the site is the marked 
site |m). The walk operator is a product of two op- 
erators, coin flip Rue ® In and the moving step S. The 
coin flip acts only on the coin space but the moving step 
S acts jointly on coin- and lattice-space as 



-) ® \x,v) 

^)®\x,y) 
i) ® \x,y) 



^)®\x + l,y), 
~^)®\x- l,y), 
i)<»\x,y + l), 
J)^\x,y-1). 



(5) 



As S involves movement only between neighboring sites, 
|a;) and \y) — > |j/± 1), it can be implemented in 

one time step. Hence Uw can be implemented in 2 time 
steps, one for = S* • {Rue In) and another for Ruc.m- 
AKR have shown that their algorithm is an instance 
of the abstract search algorithm. The operator J7vf is 
equivalent to W ■ Rue,m upto a sign, making |wc)|to) the 
effective target state \t). The walk operator W satisfies 
the required properties for the abstract search algorithm, 
within a particular subspace that is preserved by Uw- It 
is easy to check that |<i>o) is an eigenvector of W with 
eigenvalue 1. The other eigenvectors of W are 



Since the eigenstates \v~^)\xp)\Xq) are orthogonal to Tio, 
they do not matter for the algorithm and do not con- 
tribute to Ak- For even y/N, there are two eigenstates 
l*Vlv/2 %/lv/2^ of W having eigenvalue -1 within Ho- 
Since the projection of the target state on these eigen- 
states is 0{1/Vn), their contribution to Ak is negligible. 
Putting above values in Eqns. ^S^, we get 



a = e{l/VNlogN), 
|($o|a 



|(uc, m\a+)\ = 6 



y/logN 



(11) 
(12) 

(13) 



Hence, we have |$o) = \a~) + \ e), with ||e|| = 6(1/ log A''). 
After [7r/2a] = 0{^/N log N) quantum walk steps, the 
state becomes la"*") -I- |e'), with ||e'|| = 6(1/ log A). Since 
Ka+liic, m)| is Q{l/y/\ogN) and \{€'\uc,m)\ is of lower 
order 0(1/ log A^), the overlap of the final state with 
\uc,rn) is Q{l/^/\ogN). Thus we can get the \uc,m) 
state, or the |rn) state, using y/TogW rounds of quantum 
amplitude amplification. The total number of time steps 
becomes 0{y/N log N x VJogN) = 0{VN log N). 

In the next section, we show that by controlling the 
quantum step using an ancilla qubit, the coefficients apq, 
ao and Ak can be manipulated in such a way that no 
rounds of quantum amplitude amplification are required 
and the |m) state can be obtained in 0{^/N log N) time 
steps. 



I^P,) ^ Kg)\xp)\Xq) , p,gG{0,...,VAr-l}, (6) 
where \xp) - ^ E^^o"' and \x,) - 



'^^^y=o ^ e.^^'^'^'^ly) form the Fourier basis. For each 
p and q, there are four eigenvalues 1,-1 and e^**'"' with 



If 2'Kp 2'!Tq 
cosOna = - cos^= + cos— ;= 

2 V Viv Vn 



PI 



(7) 



and 



corresponding to four different vectors 
\v^^ of the coin space. W satisfies the conditions of 
the abstract search algorithm within the subspace Tio 
spanned by the eigenstates = \v^q)\Xp)\Xq) , {P,q) ^ 
(0, 0) and |$o), and |$o) is a unique eigenstate with eigen- 
value 1 within this subspace. AKR have shown that the 
operator Uw preserves this subspace. 

As shown by AKR, the vectors \vp^) are such that 

Qpq = {^pg\uc,m) = 1/V2N. We also have ao = 

{^o\uc,m) = {uc\uc){uN\m) — 1/Vn. Using these val- 
ues and Eq. ([7]) for 9pq, AKR have shown that the sums 
in Eqns. are 



e (log TV), 



^p.q l — eos6. 

i:p,«<cot2^=.e(iogiv). 



(8) 

(9) 
(10) 



III. CONTROLLED QUANTUM WALK 
ALGORITHM 

The new algorithm attaches an ancilla qubit \b) to the 
system, which controls the operations in the joint Hilbert 
space. The algorithm works in the 8A^-dimensional 
Hilbert space H — Hb He (i^ Hn, where Tib is the two- 
dimensional Hilbert space of the ancilla qubit. We use 
the subscripts b and J respectively for denoting the states 
or operations within the ancilla qubit space Tib and the 
joint Hilbert space Tij = Tie^Ti-N- (Note that Hj is the 
working space of AKR's algorithm.) 

The circuit diagram of the algorithm is shown in Fig. 



1. The initial state is |$ol 



cqw 



= 1 



un) = |1) 



\^o)akr, and it can be prepared in 0{\/N) time steps. 
The controlled quantum walk algorithm then iteratively 
applies the operator 



Uc = {Z)b ■ ciW 



{Xl)b 



ClRue.m ■ {Xs)b 



(14) 



to |$o)cqui- Note that in the figure, the operations are 
performed sequentially from left to right, while in equa- 
tions they are performed from right to left. Xg and Z 
are the single qubit gates given by 



Xs = 



cosS sinS 
—sinS cosS 



Z = 



-1 
1 



(15) 
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FIG. 1: Circuit diagram for the controUed quantum walk 
search algorithm. The Reflect and Walk boxes denote the re- 
flection operator Ruc,m and the walk operator W as deflned 
in the text. 

Let the mutually orthogonal qubit states be 

\5o) = XI\Q) = cos5\Q) + sin5\l) , 
\5i) = Xl\l) = -sin5\Q) + cos5\l) . (16) 

The operator ciRuc,m = Isn — 2|1, Uc, m) (1, u^, m\ is the 
negative of the reflection about |l)|uc)|'Tt) state. It is 
implemented by applying Ruc,m in the joint space iff the 
ancilla qubit is in |1) state. As Ruc,m can be implemented 
in one step, ciRuc.m also takes one step. Similarly, the 
controlled walk operator ciW performs a quantum walk 
W in the joint space iff the ancilla qubit is in |1) state. 
Thus Uc takes 2 time steps for implementation, one for 
ciRuc,m and another for ciW. 

For 6 = 0, the ancilla qubit is redundant and the 
new algorithm reduces to AKR's algorithm. The op- 
timal algorithm is obtained if we choose S such that 
cosS = Q{^/l/\ogN). Then measurement of the lattice 
state, after 0(v^/VToglV) iterations of the operator Uc, 
gives the desired state |m) with constant probability. The 
total complexity of algorithm is therefore 0{y/NlogN). 

To analyse the algorithm, we first show that the con- 
trolled quantum walk algorithm is an instance of the ab- 
stract search algorithm. We have {xl)h-ciRuc,m'{Xs)b — 
csiRuc,m where csiRuc,m apphes Ruc,m in the joint space 
iff the ancilla qubit is in \6i) state. Eq. ([M]) then implies 
that 

Uc = C- CsiRucrn , C = iZ)b ■ CiW. (17) 

Since csiRuc^m = Isn - 2\Si,Uc,m) {Si,Uc,m\ is equiv- 
alent to R5i,uc,m, = 2\6i,Uc,m){6i,Uc,m\ - Isn upto a 
sign, the effective target state of the algorithm is 

\ts) = \6i)K)\m). (18) 

We need to find the eigenspectrum of C and the expan- 
sion coefficients of \ts) in its eigenbasis. If |$)j is an 
eigenvector of the walk operator W with eigenvalue e*^ 
then it is easy to check that |l)6|$)j is an eigenstate of 
the operator C with the same eigenvalue e'^. Explicitly, 

e^'\l)b\'S>)j ^ e^'\l)b\^)j- (19) 



Similarly, |0)f,|$)j is an eigenvector of C with the eigen- 
value — 1 due to the Z operator. Hence the subspace 
spanned by the states |0)f,|V')j is the — 1-eigenspace of C. 
Having determined the eigenspectrum of C in terms of 
that of W, we can easily infer that C satisfies the required 
conditions for the abstract search algorithm within the 
subspace Hb ® "Hq. Moreover, this subspace is preserved 
by the operator Uc- 

To quantify the dynamics of the algorithm, we now 
calculate the quantities given by Eqns. ([TJ3]). Let ao{S), 
apq{5) and ak{6) denote the expansion coefficients of the 
target state \ts) in the eigenbasis of C . We have 

apq{5) = {6i,Uc,m\l,(^pq) = apqCosS. (20) 

where the expansion coefficients of \uc, m) in the 

eigenbasis of W , discussed in the previous section. Simi- 
larly, we get ao(^) = aocosS. Apart from these, the pro- 
jection Ak of \ts) on the —1-eigenspace of C is non-zero. 
It corresponds to the ancilla qubit being in |0) state, so 
Ak{S) = \{Si\0)\ = \sinS\. This projection was not sig- 
nificant in AKR's algorithm, but it is crucial for the new 
algorithm. 

Using these values, and Eqns. ((8 HT0|) for the sums 
occuring in Eqns. HH!]), we find that the two relevant 
eigenvectors of Uc are | ±0:5) with the eigenvalues e**"^, 
with 

as = el ^ = . (21) 

{logN+'-^)J 

The overlap of the initial state |$o)cgu) with \aj) is 

\{aJ\^o)c,^\ > 1 - e (^^^) - e{NattanH). (22) 

After T = [4^! iterations of Uc, we get the state 
Its overlap with the \Si,Uc,m) state is 

|(.r,..H«|>H-n(e(-^),l). (23) 

W e consid er the special case when cosS = 
Q{y/l/ logN). In this case, the |a^) state has a con- 
stant overlap with the desired |m) state, and hence mea- 
suring the state will give \m) with constant probabil- 
ity. Using tan'^S = e(logiV) in Eq. we find that 
as = e(l///VToglV). Putting it in Eq. we get 
|(a_^|<i>o)cq-!«| = 1 — 0(l/log^A^) so the initial state is 
very close to \aj). The required number of iterations to 
get the state \aj) is 0{1/ as) = 0{y/NlogN). Thus the 
time complexity of the algorithm is 0{^/N log N). 

If we choose cosS ^ ^/l/\ogN, then using Eq. ([23]) . 
we find that the |a^) state has still a constant overlap 
with the desired \m) state, but tan? 5 ^ logA^ and the 
number of iterations required to get the \a's) state is 
much higher than 0{-\/N \ogN). If we choose cos5 ^ 
log N , then the number of iterations required to get 
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|a^) state remains 0{y/N log A''), but this state is no 
longer close to the desired state \m) and quantum ampli- 
tude amplification is needed to get to the desired state. 
The balance is achieved when cosS = 8(-\/l/ log A^). 

IV. DISCUSSION 

We have presented a modification of the discrete time 
quantum walk search algorithm by Ambainis, Kempe and 
Rivosh for the problem of two-dimensional spatial search. 
The new algorithm solves the problem in 0{y/N log N) 
time steps and improves on AKR's algorithm by a factor 
of 0(\/log N). It can be easily generalized to the contin- 
uous time quantum walk algorithm by Childs and Gold- 
stone Q. In the continuous walk algorithm, the system 
is evolved under a time-independent Hamiltonian and 
the restriction on the Hamiltonian is that it should cou- 
ple only neighboring sites. To apply the new algorithm, 
we just attach an ancilla qubit to the Hilbert space and 
then evolve the whole system under a suitably controlled 
Hamiltonian. 

It is an open question that whether the performance of 
algorithm can be further improved. As the problem has 
a lower bound of ri{\/N) time steps it will be inter- 
esting to get an algorithm which can solve the problem 
in 0{^/N) time steps or to show that no further improve- 
ment over 0{y^N log N) complexity is possible. Within 
the framework considered here, probably 0{y/N log N) 
complexity is the best that can be achieved. The mini- 
mum eigenvalue gap of the walk operator or the Hamil- 
tonian is 0{l/y/Nlog N), so the adiabaticity condition 
demands a minimum evolution time 0{y/ N log N). Even 
the algorithms of AKR and CG evolve the system for 
0{^/N log N) time, but their final statea are not close to 
the desired state. In the new algorithm, we have intro- 
duced extra eigenstates of the walk operator by attaching 
an ancilla qubit. These extra eigenstates allow interfer- 
ence in such a way that the final state gets close to the 
desired state. 

The algorithm presented in this paper assumes a 
unique marked item, but it can be easily generalised to 
the case of multiple marked items with 0(log7V) over- 
head in computational complexity [3], making the total 
complexity of the algorithm 0(v^log'^^^ A^). In their 
paper, AKR have extended their algorithm to the case of 
two marked items (see section 6.5 of [6]), and they have 
shown that the algorithm succeeds in only 0{\/N log N) 
time steps for this case. The same extension applies 
to the new algorithm which solves the same problem 
in 0{y/N log N) time steps. Similarly, the extension of 
AKR's algorithm to the case of two-dimensional coin- 
space (see Theorem 3 of @) also applies to the new al- 
gorithm. 

Finally, we point out that the new algorithm can be 
applied to any instance of the abstract search algorithm, 
but the improvement factor may not be significant. In 
the case of higher-than-two dimensional spatial search. 



AKR's algorithm solves the problem in c^/N time steps 
where c is a constant (see Theorem 4 of Q). By using the 
new algorithm, we can improve the performance only by 
a constant factor. It can be shown that if c 1, then the 
performance can be improved by a factor of ^/c, making 
the total complexity VciV (see section III.B of [10]). For 
c — 0(1), there is not much improvement, obviously be- 
cause fl{^/N) is the lower bound on any quantum search 
algorithm. 

Note added: After this work was completed. Prof. 
Apoorva Patel pointed out to me that similar improve- 
ment in algorithm complexity can be obtained using the 
Dirac equation with a mass term [ll[ . A non-zero value 
for the mass eliminates the infrared divergence, and pro- 
vides the best performance when scaled appropriately 
with the lattice size. 

Acknowledgments: I thank Prof. Apoorva Patel for 
going through the manuscript and for helpful comments 
and discussions. 



APPENDIX: ABSTRACT SEARCH ALGORITHM 

Here, we present the analysis of the abstract search al- 
gorithm, which iterates the operator Ua — U ■ Rt on the 
state l^o) that is a unique eigenstate of U with eigen- 
value 1. Here Rt is the reflection operator about the 
target state \t) and U is required to be a real operator. 
The analysis closely follows that of AKR (see section 7 of 
6]) with the difference that we have considered the pos- 
sibility that U may have an eigenspace with eigenvalue 
— 1, referred to as the — 1-eigenspace here. We will find 
the relevant features of the eigenspectrum oi Ua, which 
are completely determined by the eigenspectrum of U 
and the expansion coefficients of \t) in the eigenbasis of 
U. As discussed in section II, the target state \t) can be 
expanded in the eigenbasis of U as 

j k 

For convenience, we use the notations ai, and 9i 
{I e {0,j,k}), for denoting the expansion coefficients 
ao,aj,ak, the eigenvectors |$o), \^k), and the eige- 
nangles Oq = 0, 6j ^ {0, tt}, 9k = tt respectively. 

We define for real A, the unnormalised vector \w\), 
whose expansion coefficients in the eigenbasis of U are 
given by {^i\wx) = aiFxiOi), FxiOi) = cot{^). We 
state some relations satisfied by function F\, which we 
will use later. These relations can be derived easily as is 
done in [fi]. 

e»('(_l + ^Fy,{e)) = e'^(l + iFx{e)), (24) 
F,(0)+F,(-^?) = ,^^fH^, (25) 
Fa(0) = cot| ; Fx(t:) = -tan|. (26) 
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As shown by AKR, if |wa) is orthogonal to \t) then the 
unnormaUsed vector |A) = \t) +i\w\) is an eigenvector of 
the operator Ua = U ■ Rt with the eigenvalue e*^. It is 
because of the special properties of the function Fx{9i). 
To see this, we note that the expansion coefficients of |A) 
in the eigenbasis of U are 

($i|A) = {^i\t) + i{Hwx) = ai{l + iFxiei)). (27) 

We have Rt\X) = —\t) + i\w\) as Rt does not alter |w^a), 
orthogonal to \t) by assumption. Hence we have 

{^i\Rt\\) = -{^i\t)+i{^i\wx) = ai{-l + iFx{ei)). (28) 

Since |$;) are eigenvectors of U, we have {^i\U • -Rt|A) = 



aie'^'i-l + iFx{9i)). Using Eqns. 



(24,27), we find it to be equal to ($(|C/a|A) = e*^ai(l + 
iFx{6i)) = e^^{^i\\). As this holds for all the basis vec- 
tors 1$;), we find that |A) is an eigenvector of Ua = U- Rt 
with eigenvalue e*"^, iff \w\) is orthogonal to \f). This 
condition is equivalent to J^i^f^xi^l) — 0- Expanding 
this sum for / = 0, j and k, and using Eq. (25) for the 
term F\{Oj) + Fx{—Oj) occuring in the sum, we find the 
condition to be 



normalized state. To show that the initial state |<I>o) is 
spanned by the eigenvectors | ± a) , we find the overlap 
of |$o) with the vector \a^). The expansion coefficients 
of the vector \a~) in the eigenbasis of U are given by 

\{^i\aZ)\ = - {^i\w-c,) = ai{F^{ei) - F.^iOi)). 

(32) 

We have |($o|a-)| = Putting / = in above 

equation, we find |($o|a~)| = ao(-Fa(0) — F_(j(0)) = 
2aocot^, and hence we need to bound ||q!~|| to bound 
|($o|q:")|. Now 



ll"«ll = jEl«'(^«W-^-a(^0)P- (33) 



^^eot| zaj ^ ^^^tanq ^To+Tj+Tk, we get |($o|a 

° sinX ^ cosX - cosOj sinX ' Zi+Zk. More exphcitly 



■hr, 



In the summation over Z, the term Tq corresponding 
to ; = is equal to Tq = \{^o\aZ)? = ^alcot^^ = 
Q{a^/ a^). Similarly the term corresponding to Z G fc 
is equal to AAltan"^^ = Q{Ala^). The term Tj corre- 
sponding to I £ j was calculated by AKR and found to 
be Tj = 6(a^ J2j /(I ~ cosOjY). Moreover, in the case 
of spatial search, tlicy have shown that Tj and Tk are 
small compared to Tq for large N. Hence, using ||q;~|| = 

= v^/ll««ll = 1- 



where Ak = v^Cfc^ ^^'^ projection of \f) state on — 1- 
eigcnspacc of U . It is easy to check that if above equation 
is satisfied for A then it is also satisfied for —A and vice 
versa. 

Let Omin be the smallest of 6j . Then as shown by AKR, 
the above equation has exactly two solutions, X = a and 
A = —a, such that |a| < O^nml'^- Moreover, the eigen- 
vectors corresponding to these eigenvalues arc relevant as 
|$o) is almost completely spanned by them, and hence it- 
eration of Ua on |$o) can be analysed by considering only 
these eigenvectors. Typically Omin is very small (in the 
case of two-dimensional spatial search, it is 0{l/^fN), 
and therefore a is very small. Writing above equation 
upto first order in a, we get 



■v2 



^ cosa — cosOj 
3 ■' 



+ 



Al 



(30) 



As shown by AKR, the first term on R.H.S. is 

2 

©(Ej T^^X which leads to 



a = e 



ao 



■COS$j 



42 



(31) 



|($o|Q!")| > 1-© a 



4 ^ 



^ (1 - cos0jy 



-e 



AW 



(34) 

Thus the state |$o) is very close to |a ) = c{\a) — \ — 
a)), where c is the normalization factor. As | ± a) are 
the eigenvectors of Ua with eigenvalues e**", we have 
([/A)«|a-) = c(e'«"|Q!) - e-'«"| - a)). After T = [7r/2a] 
iterations of Ua, we come very close to the state = 
c(|a) + |-a)). 

The last part of analysis is to calculate the overlap 
between \t) and |a+) states. Let |a+) = \a) + \ — a) be 
an unnormalized state. We have |a+) = |Q:+)/||a+|| and 

hence |(i|a+)| = As |a+) = 2|t)+i(|«;„) + |t«_„)) 

and \w±a) are orthogonal to \t), we find |(t|a+)| to be 
equal to 2. Similarly, 

11^ = l|2|i) + i{\wa) + \w-a))\f = 4 + \\wc, + w;_af . 

The expansion coefficients of the vector \wa+W-a) in the 
eigenbasis of U are given by {^i\wa+W-a) = ai{Fa{Oi) + 
F-a{0i)), and hence 



(36) 



Let I ± a) = \t) + i\'w±a) be the unnormalised eigen- 
vectors of Ua corresponding to the eigenvalues e^'". Let 

\a~) = \a) — I — a) = i{\wa) — \w^a)) be an imnormalised 
state and let \a^) = \oiu) /Wo^uW be the corresponding 



For I e {0, k}, the term Fa{Oi) + F-a{di) vanishes as 61 
is either or tt for such Vs and Fa{Oi) = —F-a{6i) for 
Oi G {0,7r}. So, all the non-vanishing terms in above 
sum correspond to I € j. This sum has been computed 
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by AKR and shown to be ©(X^j '^jCot^^). Putting it in AKR's spatial search algorithm), we get 
Eq. (35), we get 

/ / .\\"'/' \{t\a+)\ = e^m^n^ ^ l]] . (38) 

(l + e(^a>t^|jj . (37) V VV^.-.^^ // 

This completes the analysis of the abstract search algo- 
Whcn the sum is large compared to 1 (as in the case of rithm. 



[1] L.K.Grover, Phys. Rev. Lett. 79, 325 (1997). 

[2] P.Benioff, AMS Contemporary Mathematics, 305, pp. 
1-12 (2002). e-print quant-ph/0003006 

[3] S. Aaronson and A. Ambainis, Proc. 44th IEEE Sympo- 
sium on Foundations of Computer Science, pp. 200-209 

(2003) . e-print quant-ph/0303041 

[4] N. Shenvi, J. Kempe, and K.B. Whaley, Phys. Rev. A 

67, 052307 (2003). 
[5] A. Childs and J. Goldstone, Phys. Rev. A 70, 022314 

(2004) . 

[6] A. Ambainis, J. Kempe, and A. Rivosh, Proc. 16th 
ACM-SIAM SODA, pp. 1099-1108 (2005). e-print 
|quant-ph/0402107| 



[7] A. Childs and J. Goldstone, Phys. Rev. A 70, 042312 
(2004). 

[8] C. Bennett, E. Bernstein, G. Brassard, and U. Vazirani, 
SIAM Journal on Computing, 26, 1510 (1997). e-print 
quant-ph/9701001 
[9] G. Brassard, P. Hoyer, M. Mosca, and A. Tapp, AMS 
Contemporary Mathematics, 305, pp. 53-74 (2002). e- 
print quant-ph/0005055 

[10] A. Tulsi, A general framework for quantum search algo- 
rithms, manuscript in preparation. 

[11] A. Patel, Poster presented at QIP-2008, New Delhi, De- 
cember 2007, unpublished. 



